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Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress 
tensor and its parameterization in terms of turbulent viscosity and A-effect. The turbulence is forced in an anisotropic 
fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis 
inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified 
convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, 
such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here. 
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1 Introduction 

The Reynolds stress, described by the correlation of fluc- 
tuating velocity components, Qij — UiUj, is one of the 
most important generators of differential rotation in stars 
(Riidiger 1989). These stresses have been studied with the 
help of 3D convection simulations (e.g. PuUckinen et al. 
1993; Chan 2001; Kapyla et al. 2004; Riidiger et al. 2005). 
These results have revealed some surprising features such 
as the peaking of the horizontal stress Q^y very close to 
the equator, and a positive (outward) flux for rapid rotation. 
Both of these results are at odds with theoretical considera- 
tions (Kitchatinov & Riidiger 1993). Furthermore, disentan- 
gling of the diffusive (turbulent viscosity) and non-diffusive 
(A-effect) parts of the stress is difficult from convection sim- 
ulations. 

Here, we present preliminary results from anisotropic 
homogeneous, isothermal, non-stratified turbulence simu- 
lations in which diffusive and non-diffusive effects can be 
studied separately. Imposing a linear shear flow on top of 
isotropically driven turbulence allows the study of turbulent 
viscosity without A-effect. On the other hand, using a spe- 
cial form of forcing, anisotropic homogeneous turbulence 
can be generated. Rotation is added to study the A-effect. A 
simple analytical closure model, based on the minimal tau- 
approximation (hereafter MTA, see e.g. Blackman & Field 
2002; Brandenburg et al. 2004), is used to compare with 
simulations in the cases with rotation. 

2 The models 
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-<V Inp - 2 X M + /force + /vise , (2) 



using an isothermal equation of state characterized by sound 
speed Cg in a fully periodic cube of volume (27r)'^. Here, 
D/Dt = d/dt — M • V is the advective derivative, p is the 
density, m is the velocity, and ft = ilo{— smd,0,cos9)^ 
is the rotation vector By virtue of the periodic boundaries 
mass is conserved and the volume averaged density has a 
constant value of p = po- In the present study we use an 
anisotropic forcing function in Fourier space according to 



(force) 



= (/o% + fl COS^efc Z,Zj)f^ 



(iso) 



(3) 



where /o is the amplitude of the isotropic part and /i the 
anisotropic one, 9fc is the angle between the vertical direc- 
tion and the wave vector k, and z is the unit vector in the 
vertical direction. Details of the isotropic part of the forcing 
are given, e.g., in Brandenburg et al. (2004). 
The viscous force is given by 
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where v is the viscosity and 
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(5) 



is the rate of strain tensor The simulations were made with 
the Pencil-Cod^H 

In the minimal tau-approximation one solves for the 
time derivative of a quantity instead of the quantity itself, 
i.e. in the present case 



In the 3D simulations we solve the set of equations 
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Qij iliUj 



(6) 



( 1 ) where dots denote time derivatives and the overbars volume 
averages. Inserting the Navier-Stokes equations into Eq. ^ 



http : / /www .nordita. org/ software /pencil-code/ 



®W]LEY 

JiTiterScieTice* 



© 0000 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



P.J. Kapyla & A. Brandenburg: Turbulent viscosity and A-effect from numerical turbulence models. 



2 







Turbulent 


viscosi 


ty 


100 




^ ■ 










slope = 


1.52 ... 

A' 


256'' 


\ 






128'' 






i" 


.4' 
64'' 




1.52 Re 


10 











10 100 
Re 

Fig. 1 The ratio of turbulent to molecular viscosity, ft /i^, 
as a function of the Reynolds number. The resolution used is 
denoted alongside each point. The error bars are estimated 
from a modified mean error of the mean; see, e.g. Eq. (30) 
of Kapyla et al. (2004). 

and assuming a high Reynolds number so that the viscous 
terms can be neglected, one arrives at 

Qij = —2ejki^kQii — ^Eiki^kQji + 

Uifj + Ujfi + Tij , (7) 

where Tij denotes the triple correlations. In the present case 
we assume that the pressure terms are subsumed in the 
triple correlations, which is reasonable for calculating the 
A-effect, but not valid for calculating turbulent viscosity. 
The basic assumption of MTA is that the triple correlations 
can be presented in terms of the quadratic ones 

Tij = -T~^Qij , (8) 

where r is a relaxation time. In the statistically steady state 
without rotation we have 

Qy =t(^ + ^) , (9) 

which allows us to express the forcing in terms of the 
Reynolds tensor for the non-rotating case (denoted by the 
superscript zero). We employ the same amplitudes as those 
found in the 3D simulations. The only free parameter in the 
MTA-model is then the Strouhal number 

St = TUrmskf , (10) 

where u^-^s is the rms-value of the fluctuating component of 
the velocity and fcf « 5 is the mean forcing wave number 
In the present study we use St = 2 which reproduces the 
same trend as a function of rotation as the 3D simulations 
(see Fig.|2]i. 

3 Results 

3.1 Turbulent viscosity 

In order to study the turbulent viscosity, a large scale linear 
shear flow U = (0, Sx, 0) is imposed upon the system. In 



this system the homogeneity of the turbulence is preserved 
by using the shearing box approximation (e.g. Hawley et al. 
1995). The turbulent viscosity can now be computed from 

i^t = -2Q„j/S. (11) 

where S = —0.1 is used in the present study. Fig. [T] shows 
z/t /j^ as a function of Reynolds number. 

Re = ^ . (12) 

The absolute value of the turbulent viscosity stays almost 
constant as a function of Reynolds number and thus the ra- 
tio lyt/iy increases almost linearly as a function of Re except 
for the largest Reynolds number run. Higher resolution sim- 
ulations are needed in order to clarify the behaviour in the 
high Reynolds number regime. 

3.2 A-effect 

In anisotropic turbulence under the influence of rotation 
non-zero off -diagonal Reynolds stresses are generated ac- 
cording to (Riidiger 1989) 

Qij = Aijkflk + diffusive terms. (13) 

Considering the computational domain as a small rectangu- 
lar part of a sphere, we note that Qxy describes horizontal 
(latitudinal), and Qy^ vertical (radial) transport of angular 
momentum. In the present study the vertical {z) direction is 
taken to be the preferred one. Furthermore, we use /i = 0.2 
and fo = which results in turbulence that is dominated 
by the z-component which, bearing in mind the geometry 
of the system, is also the case for convection. The ratio of 
vertical to horizontal turbulence intensity is w^/w^ « 2.2, 
where = (u^ + Uy)/2. A resolution of 32^ was used in 
the present calculations with Re « 7. The Reynolds num- 
ber dependence of the A-effect is weak if Re is sufficiently 
large (not shown). This issue will be studied further in a fu- 
ture publication. 

The left hand panels of Fig. |2] show the results from di- 
rect 3D calculations. The rotational influence is quantified 
by the Coriolis number 

Co = 2r!oKmsfcf)"^ (14) 

Note that in comparison to the commonly used definition 
in convection simulations, our Coriolis numbers are smaller 
by a factor of 2tt. The main feature of the horizontal stress 
is that it is always positive and peaks around 30 degrees 
latitude for all Coriolis numbers studied so far This result 
is in stark contrast to the convection calculations where Qxy 
peaks always very near the equator (e.g. Chan 2001 ; Kapyla 
et al. 2004; Hupfer et al. 2005). 

The vertical stress is predominantly negative with a 
maximum at the equator Although there seems to be a 
regime in which the vertical stress is positive for interme- 
diate rotation rates, the sign remains negative for the most 
rapid rotation cases studied so far This feature is also at 
odds with convection simulations, which exhibit positive 



© 0000 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumal.or^ 



Astron. Nachr. / AN (0000) 



3 



Qxy/uLs(3D calculation) 




0.00 
-0.02 
-0.04 
-0.06 



QyyuLs(3D calculation) 













" I, 


- Q = 


0. 02, Co = 0.05 


Q = 


0.1, Co = 0.3 \ 


n = 


0.5, Co=1.5 ' 


■ n = 


2.0, Co = 5.3 







30 



60 



90 



Qxy/uLs(MTA-closure) 




Qyz/uLs(MTA-closure) 



0.00 
-0.02 
-0.04 
-0.06 





.\ ^^-^^^ 




\ — - — 




\ 'v, ^ - 

- — n '= 0.02 ' 




n ='0.1 




1 n = o'.5. 




n = 2.0 





30 



60 



90 



Fig. 2 The non-diffusive Reynolds stresses (i.e. A-effect) Q^y (upper panels) and Qyz (lower panels) from direct 3D 
simulations (left panels) and the MTA-closure (right panels) as functions of latitude and rotation rate. The error bars are 
defined similarly as in Fig.[T] 



Qyz for rapid enough rotation (Kapyla et al. 2004; Chan 
2007, private communication). 

The MTA-model results are shown in the right panels of 
Fig- El Qualitatively the results match the numerical simu- 
lations rather well. The sign and latitude distribution of Q^y 
is reproduced fairly well. The amplitude, however, is clearly 
too large. Similar conclusions can be drawn from the results 
for Qyz, although there the latitude distribution from the 
closure model tends to show a persistent maximum at mid- 
latitudes which is not observed in the 3D simulations. Also 
the amplitude is too large by almost a factor of two. 

A possible explanation of the discrepancies between the 
numerical results and the closure model is that the latter 
does not take isotropizing effects properly into account. 
This can be due to an improper treatment of the pressure 
terms in the closure model. 

4 Conclusions 

Shear flow turbulence simulations show that the ratio of tur- 
bulent to molecular viscosity increases linearly up to Re w 
30 with I't/i' ~ 1.5 Re. For the largest Reynolds number 
the scaling seems somewhat shallower but the present data 
is not yet sufficient to substantiate this. The A-effect from 
homogeneous, anisotropic turbulence does not exhibit the 
puzzling features found in convection simulations. Further 



study is required in order to understand which of the ne- 
glected physics is responsible for the lack of these features. 

The MTA-closure is able to reproduce many of the qual- 
itative aspects of the simulation results including a maxi- 
mum of the horizontal stress at about 30° latitude, with its 
largest value for Co « 0.5. In the model, the vertical stress 
can have a maximum away from the equator for Co <; 0.2, 
which is not seen in the simulations. Nevertheless, both 
simulations and model have the largest vertical stress for 
Co w 0.3. However, the model generally overestimates the 
magnitudes of the stresses. More detailed analysis of the 
simulation and closure results will be presented in a future 
publication. 
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